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ABSTRACT 

o. 

£N) _ Context. C3 is the smallest pure carbon chain detected in the dense environment of star forming regions, although diatomic C2 is 
detected in diffuse clouds. Measurement of the abundance of C3 and the chemistry of its formation in dense star forming regions 
has remained relatively unexplored. 

Aims. We aim to identify the primary C3 formation routes in dense star forming regions following a chemical network producing 
species like CCH and C-C3H2 in the star forming cores associated with DR21(OH), a high mass star forming region. 
Methods. We have observed velocity resolved spectra of four ro-vibrational far-infrared transitions of C3 between the vibrational 
ground state and the low-energy V2 bending mode at frequencies between 1654-1897 GHz using HIFI on board Herschel, in 
DR21(OH). Several transitions of CCH and C-C3H2 have also been observed with HIFI and the IRAM 30m telescope. Rotational 
' temperatures and column densities for all chemical species were estimated. A gas and grain warm-up model was used to obtain 
. estimates of densities and temperatures of the envelope. The chemical network in the model has been used to identify the primary 
\^ • C3 forming reactions in DR21(OH). 

Results. We have detected C3 in absorption in four far-infrared transitions, P(4), P(10), Q(2) and Q(4). The continuum sources 
MM1 and MM2 in DR21(OH) though spatially unresolved, are sufficiently separated in velocity to be identified in the C3 spectra. 
All C3 transitions are detected from the embedded source MM2 and the surrounding envelope, whereas only Q(4) & P(4) are 
q ■ detected toward the hot core MM1. The abundance of C3 in the envelope and MM2 is ~ 6 x 1CP 10 and ~ 3 x 10~ 9 respectively. 
Uh " For CCH and C-C3H2 we only detect emission from the envelope and MM1. The observed CCH, C3 and C-C3H2 abundances are 
1 most consistent with a chemical model with tih 2 ~ 5xl0 6 cm -3 , post-warm-up dust temperature, T max =30 K and a time of ~ 
. 0.7-3 Myr. 

Conclusions. Post warm-up gas phase chemistry of CH4 released from the grain at t ~ 0.2 Myr and lasting for 1 Myr can explain the 
observed C3 abundance in the envelope of DR21(OH) and no mechanism involving photodestruction of PAH molecules is required. 
The chemistry in the envelope is similar to the warm carbon chain chemistry (WCCC) found in lukewarm corinos. We interpret 
^ , the observed lower C3 abundance in MM1 as compared to MM2 and the envelope to be due to the destruction of C3 in the more 
vo. evolved MM1. The timescale for the chemistry derived for the envelope is consistent with the dynamical timescale of 2 Myr derived 
. for DR21(OH) in other studies. 

00 • 
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1. Introduction ward the Sgr B2 star forming region and IRC +102 16 by 

ICernicharo et all (2000). G iesen et al.l (|200lh discussed new 
Small carbon chains are important for the chemistry of laboratory data on the vibrational spectrum of C 3 in its 
stellar and interstellar environments for several reasons: i ow _f reque ncy bending mode and re-visited the first iden- 
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they are ubiquitou s throughout the interstellar medium tific ation of the R(2) l ine in absorption toward Sgr 

C3| flAdamkovics et al.| [2003|), they arc likely to participate B2 (|Van Orden et all 119951 ). The abundance and excita- 

in the formation of long carbon chain molecules, and tion of C 3 with a large range of rotational temperatures 

they are products of photo-fragmentatio n cascades o f poly- i n tra nslucent clouds have also been determined convinc- 

cyclic aromatic hydrocarbons (PAHs) (|Radi et al.| |l988j ingfy ([Maier et alj[200lllRoueff et al.ll2003:IOka et alj2003t 

l Petv et aL l ™)- Tri atomic carbon, C 3 , was first ten- | A damkovics et~all [200l through observations at optical 

tatiyely identified m interstellar gas by | Van Orden et al. | wavfi i fiTlrt > ia . lOalWdiuov et al l Hw(Y% also demonstrated 



tatiyely identified m mterstehargas by I Van CJrden et al. | wave len g ths. lGala Z utdinov et atl $2(M) also demonstrated 
H995D and| Haffner fc Meyer| ([1995|) . The mid-infrared spec- fh . , hp Cq RhllllHRT1 „ P 1s rfilatfiH UP ; thpr tn ^terstellar red- 



that the C3 abundance is related neither to interstellar red- 
trum of C 3 (1* antisymmetric stretching mode) was mea- dcning nQr to the intensitics of diffuse int erstellar bands, 
sured in t he circumstel l ar en velope of CW Leo (IRC 

+10216) bv iHinkle et al.l (|l988tl . and in low-resolution in- DR21(OH) lies 2' North of the H 11 region DR21, in 
terstcllar absorption in the far-IR [y^ bending mode) to- the Cygnus X H 11 complex. The distance to DR21(OH) 
has recently been accurately determined by trigonometric 



* Herschel is an ESA space observatory with science instru- parallax measure ments of its associ ated methanol masers 
ments provided by European- led Principal Investigator consor- as 1.50t ; 07 kpc (|Rygl et alj 1201 if ). Interferometric and 
tia and with important participation from NASA. high resolution single-dish observations both in contin- 
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Table 1. Properties of the sources in DR21(OH). 



Component 


q(2000) 


£(2000) 


vlstl 
(km s" 1 ) 


Mass" T d 
(M ) 


N(H 2 ) 6 
(10 23 cm" 2 ) 


< n{H 2 ) > 
(10 7 cm" 3 ) 


References 


MM1 


20 h 39 m 0F0 


42°22'48" 


-4.1±0.3 


88 58 


4.5 


4.1 


1,2 


MM2 


20 h 39 m 00:4 


42°22'43'.'8 


-0.7±0.3 


143 30 


7.4 


6.8 


1,2 


Envelope 






-3.1±0.3 


118 c 


3.5 


0.3 


1,3,4 



a All estimated for a distance of 1.5 kpc 
b Within a 15" beam 

c Mass of the enve lope derived by removing the contributions of MM1 and MM2 from the total mass of the region estimated by 
iMotte et all 1120071). 

References: (1) iMangum et al.l <|199in . (2) iMangum etHI (j!992l ). (3) IWilson fc MauersbergeJ (| 19901 ). f4 lMotte et al.l (|2007T >. 
Table 2. Spectroscopic parameters for the C3 & CCH transitions observed with HIFI/Hcrschcl 



Species 


Transition 


Frequency 


A-coeff 




Ei 


Beam Size 






(MHz) 






(cm- 1 ) 




C 3 , (J,v) 


(9,1) - (10,0) P(10) 


1654081.66 


2.38 10" 


-a 


47.3 


13" 




(3,1) - (4,0) P(4) 


1787890.57 


2.72 10" 


-3 


8.6 


12" 




(2,1) - (2,0) Q(2) 


1890558.06 


7.51 10" 


-3 


2.6 


11" 




(4,1) - (4,0) Q(4) 


1896706.56 


7.58 10" 


-3 


8.6 


11" 


CCH, Nj, F 


6l3/2,7 ~ 5n/2,6 


523971.5704 


4.58 10" 


-4 


43 


40" 




6l3/2,6 ~ 5n/2,5 


523972.1630 


4.53 10" 


-4 


43 


40" 


CCH, N JiF 


6ll/2,6 ~ 5n/2,5 


524033.9075 


4.51 10" 


-4 


43 


40" 




6ll/2,5 _ 5g/ 2j 4 


524034.5305 


4.43 10" 


-4 


43 


40" 



* Herschel HIFI beamsizes are taken from lRoelfsema et al.l (|2012T ). 

Table 3. Spectroscopic & Observational parameters for the transitions of CCH and C-C3H2 observed with IRAM 30 m. 



Species 


Transition 


Frequency 


A-coeff 


E, 


#FWHM 


^eff 


Befl 






(MHz) 




(cm' 1 ) 


(") 






CCH, Nj, F 


13/2,1 ~ 0i/ 2j i 


87284.1050 


2.60 10~' 


0.0015 


29" 


0.95 


0.75 




13/2,2 ~~ 0i/ 2j i 


87316.8980 


1.53 10" 6 


0.0015 


29" 


0.95 


0.75 




13/2,1 ~~ 0i/ 2 ,0 


87328.5850 


1.27 10~ 6 


0.0015 


29" 


0.95 


0.75 


C-C3H2, (J Ka ,K c ) 


2l,2 - lo,i 


85338.8930 


2.32 W~ b 


1.7 


29" 


0.95 


0.75 



a F e g Forward efficiency and b B e g Beam efficiency of the telescope. 



uum and m olecular lines have shown multiple peaks in 
DR2 1 (OH) (|Wilson fc Mauersbergeilll990l: IMangum et al.l 

[T991T [199a 



Chandler et all 119931) . IMangum et al 



(|1992D 

resolved the main DR21(OH) peak into two sources, MM1 
and MM2, and MM2 into 2 sub-sources, MM2-A and MM2- 
B. Synthesis imaging of NH3 emission from the region 
has clearly resolved MM1 and MM2, having radial veloc- 
ities — 4.1±0.3 kms" 1 and — 0.7±0.3 kms -1 , respectively 
(|Mangum et al.l fl992t ) . Both the star- forming cores MM1 
& MM2 in DR21(OH) are young, having no visible H 11 
region and only very weak continuum em ission at the cen- 
timeter wavelengths (| Argon et al.l l200"ol ). Dust continuum 
observations suggest that MM1 is the brighter source (L = 
1.7 x 10 4 L Q ; ZAMS B0V star) showing evidence of star 
formation, whereas MM2 (L — 1.2 x 10 3 L Q ; early B star) 
though more massive is fainter and most likely at an earlier 
stage of evolution. The region is characterized as a high- 
mass star-forming region due to the detection of centimeter 
and millimeter maser e mission from numerous transitions 
e.g., H 2 Q, OH CH3OH fArava et al.ll2009HFish et al.l l2005t 
IMangum et"ai1ll992l:|Plambeck fc Mentenlll990l) . MM1 also 
contains ground-state OH masers at 1.6 GHz and 6.7 GHz 
class II CH3OH masers near the peak of the dust and cen- 



timeter continuum emission (|Argon et al.ll200ri iFish et al.l 
120051 : iRvgl et all l201ll ) . In a recent interferometric study 
of the continuum e mission at 1.4 mm from DR21(OH), 
IZapata et ail (|2012h have resolved MM1 and MM2 into 9 
compact sources. Five of the compact sources are associated 
with MM1 and four with MM2. Two of the compact sources 
associated wit h MM1 (SMA6 and SMA7) seem to show hot 
core activity. IWhite et all (|2010D have observed multiple 
transitions of CO from 4-3 to 13-12 using SPIRE/Herschel 
to conclude the presence of two components at 80 K and 
180 K from a rotation temperature analysis. From a LVG 
analysis these authors obtain Tki n = 125 K and nn 2 = 
7xl0 4 cm -3 . We summarize selected properties of all the 
emission components in DR21(OH) as available from lit- 
erature in Table [T] All values in Table Q] correspond to a 
distance of 1.5 kpc. 

The Heterodyne Ins trument for the Far-Infrared (HIFI; 



Ide Graauw et al. 20101) on boa rd the Herschel Space 
Observatory (jPilbratt et al.ll2010l ). with its broad frequency 
coverage, high sensitivity and spectral resolution has pro- 
vided for the first time the opportunity for a systematic 
study of carbon chain molecules such as C3 through prob- 
ing several ro-vibrational lines at full spectral resolution. 
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The results presented here are a part of the PRISMAS 
("PRobing IntcrStellar Molecules with Absorption line 
Studies) Key Program (jGerin et al.ll2010D . 

The primary aim of the paper is to understand the C3 
formation mechanism. Observations of diffuse interstellar 
gas at optical and IR wavelengths have shown a strong cor- 
relation between the column densities of C3 and C2. This 
suggests that C3 and Co are i n the same chain of chem- 
ical reactions (|Oka et al.ll2003D . The formation routes for 
C3 starting from C2 also involved the production of CCH 
and C-C3H2. However no transitions of C2 are available at 
longer wavelengths (far-infrared and longer) which can be 
used to detect C2 in den se star forming regions where C3 
is detected in absorption (jMookeriea et al.l[2010Tl . Thus, to 
understand the formation pathway of C3 we compare the 
abundances of C3, CCH, and C-C3H2, with results of chem- 
ical models of dense star forming cores. 

2. Observations 

2.1. Herschel Data: C 3 & CCH 

We have observed in DR21(OH) four lines of the v-i fun- 
damental band, P(4), Q(2), Q(4) and P(10), of triatomic 
carbon using the upper sideband of the HIFI bands 7a, 
7b and the lower sideband of band 6b of the HIFI re- 
ceiver. All C3 spectra were observed at a single position 
in DR21(OH) with coordinates a(J2000) = 20:39:01.00 
(5(J2000)= 42:22:48.0, in the dual beam switch mode. 
The observations were carried out on 2010 June 25, and 
November 25 to 27. Th e details of the ob s ervat ions and data 
reduction are given in iMookeriea et all (|2010t ). All spectra 
were smoothed to resolution between ~ 0.16 to 0.18 kms -1 
and the rms noise level for the spectra is 0.03 to 0.04 K. 
All the observed C3 lines are at frequencies for which HEB 
mixers are used in HIFI, and we have noticed a variation 
in the level of the continuum observed for different LO set- 
tings and polarizations for a particular spectral line. In the 
results presented here there is only a shift in the continuum 
level, with little or no change in the depth of absorption. 
We have verified that this shift in the continuum level does 
not arise due to differences in the pointing of the telescope 
or the notorious HEB-band standing waves. We find the un- 
certainty in the continuum level to be the most pronounced 
for the P(4) line at 1787 GHz. We have evaluated the ef- 
fect of these variations and included them in our estimate 
of uncertainties of the column densities (Sec. 3.1). 

Fig. [1] shows the observed spectra normalized to 
the (single sideband) continuum level. The observed C3 
absorption features occur around the systemic velocity 
(~ — 1 km/s) of an infal ling subfilament on DR21(OH) 
(see ISchneider et al.l l2010t ) and hence are associated with 
the source only and not with foreground material. The 
PRISMAS observations are not sensitive enough to detect 
C3 lines in absorpti on in the diffuse gas a long the line of 
sight to DR21(OH) (jMookeriea et alJl2010T >. 

As part of the PRISMAS program, in addition to the 
target lines detected in absorption towards DR21(OH), a 
large number of emission lines falling within the multi- 
ple observed bands were also detected. Here we use the 
two spectral lines of CCH at 523.97 and 524.03 GHz in 
DR21(OH) in HIFI Band la. 

Table [2] summarizes the spectroscopic parameters of the 
species and transitions observed with Herschel. 



1 1 1 1 1 1 1 1 1 1 |' 1 1 'I 1 1 1 ' 1 1 1 1 1 1 




-15 -10 -5 5 10 

Velocity (km/s) 



Fig. 1. HIFI observations of the C 3 lines in DR21(OH). 
The spectra are corrected for the DSB continuum and nor- 
malized to the SSB continuum level. The Q(2), P(4) and 
P(10) spectra are shifted consecutively by +0.5, so that 
the plots do not overlap. The red line shows the simultane- 
ous fit with 3 Gaussian velocity components (at —0.7, —2.9 
and —4.2 kms -1 ). Note that the 3rd component is only de- 
tected above the noise in Q(4) and P(4). The "negative" 
absorption bumps at this velocity fitted in Q(2) and P(10) 
is within the noise. 



2.2. I RAM 30m data: CCH & c-C 3 H 2 

The ground state transition of CCH at 87.316 GHz (N=l-0, 
J=3/2-l/2, with three hyperfine components) and the 2^2- 
I0.1 transition of C-C3H2 at 85.338 GHz were observed by 
iGerin et all (|2011( ) using the IRAM-30m telescope in 2006, 
August and December. Table [3] presents both spectroscopic 
and observational parameters of the lines observed with 
IRAM 30m. The primary aim of these observations was the 
study of CCH and 0C3H2 in the diffuse and translucent 
clouds along the line of sight to DR21(OH). However, as do 
the HIFI data, these observations also detect the emission 
lines arising from the source DR21(OH) itself. 
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20 h 39 m 02" 00' 



Fig. 2. The 3 mm PdBI continuum image of DR21(OH) in 
color and black contours overlaid with yellow contours of 
450 /jm image observed with JCMT SCUBA. Also over- 
laid on the continuum image are the FWHM beam sizes 
for spectroscopic observations. The dashed lines show the 
HIFI (12" (black) & 43" (white)) and IRAM 30m (29", 
red) beams. The white cross denotes the center of the HIFI 
beams. 



2.3. Position of beams relative to the continuum sources in 
DR21(OH) 

Figure [2] presents the 12" (C 3 ) and 43" (CCH at 522 GHz) 
beams of HIFI and the 29" IRAM 30m beams overlaid on 
the images of 3 mm and 450 /im continuum emission from 
DR21(OH). The 3 mm continuum e mission (color fc black 
contours) was observed with PdBI (|Csengeri et al.l 1201 ll ) 
and the 450^m (yellow contour) continuum data were ob- 
tained from the JCMT archives. We see that the 3 mm 
interferometric map (restored to a beam of 2") clearly re- 
solves the continuum sources MM1 and MM2 and also 
shows evidence for extended emission from the envelope. 
The 450 /im JCMT data, at a resolution of 8" shows that 
the continuum emission is primarily dominated by the sur- 
rounding extended envelope. The C3 beam of HIFI covers a 
significant part of the source MM2 as well. The 29" IRAM 
30m beam almost completely covers the larger region sam- 
pled by the 522 GHz HIFI beam. Thus Fig. [2] shows that 
the observing beams for all the spectral observations pre- 
sented in this paper cover regions including MM1, MM2 
and the surrounding envelope. The spectral indices of the 
continuum emission derived so far from the 3 mm (PdBI) 
and 1.4 mm fSMA IZapata et"aT1l2012t) data are not con- 
sistent with thermal emission from dust. It is thus difficult 
to characterize the contribution of MM1 and MM2 to the 
continuum emission. However, at the beamsizes of our C3 
observations at THz frequencies it is likely that the con- 
tinuum is dominated by the emission from the extended 
envelope and hence variation due to the contribution of the 
individual sources is likely to be small. We thus assume that 
all the C3 absorption is of the same continuum, produced 
primarily from the extended envelope. 



3. Results 

3.1. C3 column densities 

The observed C3 absorption spectra show evidence for mul- 
tiple velocity components (Fig. [T]). The <2(4) transition 
clearly shows three velocity components, while the other 
transitions distinctly show two components. In a global 
fitting scheme we have assumed each C3 transition to be 
having contributions from three velocity components and 
then fitted all the observed transitions simultaneously using 
multi-component gaussians. The parameters being fitted si- 
multaneously in each transition are the common i^lsR; Ai> 
for each of the three velocity components and the strength 
of the absorption dips, separately for each line. Thus, we 
have assumed all three velocity components to be present 
in all transitions even if the components are not explicitly 
seen, since we do not attribute the non-detection to any 
physical or chemical reason rather to the lower strength of 
the transition relative to the achieved rms of our observa- 
tions. All the four spectra show evidence for two major ve- 
locity components centered at vlsr = —2.85 ±0.04 kms -1 
and i>lsr = —0.68 ± 0.02 kms -1 , with the latter com- 
ponent being stronger. The third velocity component de- 
tected only in Q(4) and marginally in the P(4) spectra is 
at —4.2 ±0.2 kms" 1 . We have checked that for a particular 
transition in which the third component is undetected the 
results of fitting do not change significantly for the well de- 
tected velocity components by the assumption of the pres- 
ence of the third component. Table 0] presents the single 
side-band continuum levels (estimated from double side- 
band temperatures assuming a side- band gain ratio of 1) 
and the results of the Gaussian fit along with their uncer- 
tainties. 

We compare the observed velocity components at —0.7, 
—2.9 kms -1 and —4.2 kms -1 with the known velocities in 
the region. The LSR velocities of the sources MM1 and 
MM2 are —4.1 and —0.7 kms -1 , respectively, and at a res- 
olution of 12" the peak of C 18 emission is observed to have 
a dlsr of —3.2 kms -1 , which most lik ely arises from the en- 
velope surrounding MM1 and MM2 (jManeum et al.l 119921: 
IWilson fc Mauersberger|[T990[) . Our observation is centered 
on MM1; MM2, located 8"3 away from MM1, lies within 
the HIFI beam. Thus, although we can not resolve MM1 
and MM2 spatially, the velocities of the two sources and 
the envelope differ sufficiently from one another to allow 
identification of the three components in the C3 spectra. 
We detect C3 absorption due to MM2 and the envelope 
in all transitions, while MM1 appears only in the stronger 
Q(4) transition and marginally in P(4) transition. Based 
on the J = 4 column densities derived for the —4.2 kms -1 
component assuming T rot between 50-70 K we estimate the 
expected absorption depths for the J = 2 level to be only 
10% and even less for the J = 10 level and this is com- 
parable to the uncertainty in our measurements. Thus the 
non-detection of the —4.2 kms -1 component in the Q(2) 
and -P(IO) is consistent with the expected strengths of these 
lines. 

Taking into consideration the uncertainties in the mea- 
surements described above, we estimate the state-specific 
column densities to be ~ 10 14 cm -2 and ~ 10 13 cm -2 for 
the —0.7 kms -1 and —2.9 kms -1 , components, respectively 
(Table H]). For Q(4) & P(4) the additional absorption fea- 
ture at —4.2 kms -1 corresponds to a J — 4 column density 
of ~ 10 13 cm -2 . Given the marginal detection at P(4) and 
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Table 4. Single sideband continuum antenna temperatures 
(T c ) and parameters derived from simultaneous Gaussian 
fitting of all the line profiles and the lower state C3 col- 
umn densities (N\) estimated from the fitted intensities. 
The errorbars in N\ include errors due to both uncertainty 
in the continuum level and the line fitting procedure. The 
velocity components are as follows: Component 1: v ccn = 
-0.68±0.02 kms -1 , AV = 1.85±0.51 kms" 1 , Component 
2: Vccn = -2.85 ± 0.04 kms" 1 , AV = 0.86±0.08 kms" 1 
and Component 3: v con = —4.2 ± 0.2 kms -1 , AV = 
1.4±0.4 kms -1 . Velocity component 3 is detected only 
above the noise level in P(4) and Q(4). 



12.6 



Transition 


T c 
[K] 


Vccn 

[km/s] 


J rdv 
[km/s] 


Ni 

[10 14 cm- 2 ] 


P(10) 


7.1T0.4 


-0.7 
-2.9 
-4.2 


0.60T0.07 
0.17±0.04 


1.2T0.08 
0.3T0.06 


P(4) 


7.0T1.3 


-0.7 
-2.9 
-4.2 


1.00T0.06 
0.08T0.03 
0.07T0.03 


2.5T0.1 
0.2T0.07 
0.17T0.08 


0(2) 


7 1 +0.6 

' ■ J --0.5 


-0.7 
-2.9 
-4.2 


1.46T0.05 
0.19±0.03 


1.2T0.06 
0.2±0.04 


Q(4) 


6 7+ - 3 


-0.7 
-2.9 
-4.2 


1.79T0.04 
0.33±0.03 
0.16T0.03 


1.5±0.06 
0.3±0.03 
0.1±0.06 


Total a 




-0.7 
-2.9 




9.1T2.3 
2.2T0.4 



a Estimated using rotation diagram (Fig. 



T — I — I — I — I — I — f~ 



-0.7 km/s 




13.6 



I — I — I — I — 1~ 



20 40 60 
E (K) 

Fig. 3. Rotation diagram of C3 for the velocity components 
at —0.7 kms -1 and —2.9 kms -1 . The errorbars correspond 
to an estimated error of 10% in the state specific column 
densities and are conservative estimates which include the 
uncertainties due to continuum measurements as well as 
the noise in the spectra. 



Table 5. Parameters derived from fitting the CCH (ob- 
served with HIFI and IRAM) and c-C 3 H 2 lines with LTE 
models using XCLASS. Parameters without errorbars were 
held constant for the fitting. 



Species 


Size 


T 


N 


Av 


VLSR 




11 


K 


10 14 cm- 2 


kms -1 


kms -1 


CCH-87 


15 


110±15 


19.7T3.0 


1.6T0.1 


-0.6 ±0.1 




30 


9±2 


23.4±1.3 


3.5±0.1 


-2.7 ±0.1 




15 


111±10 


60.3±4.0 


1.7±0.1 


-4.2 ±0.1 


CCH-524 


30 


33±5 


1.2T0.3 


5.5±0.2 


-2.8 ±0.2 




15 


101±11 


1.4±0.3 


4.5±0.2 


-4.6 ±0.2 



non-detection at Q(2) we propose the J = 4 column den- 
sity for the —4.2 kms -1 component to be an upper limit. 
Using the J = 2. 4 and 10 levels we construct rotation di- 
agrams and obtain rotation temperatures of 76±20 K and 
48±6 K for the velocity components at —2.9 kms -1 and 
—0.7 kms -1 respectively (Fig. [3]). The total C3 column den- 
sities for the —2.9 kms -1 and —0.7 kms -1 components are 
(2.2±0.4)xl0 14 cm -2 and (9.1±2.3) x 10 14 cm" 2 , respec- 
tively. Thus, the —0.7 kms -1 component corresponds to 
the colder and higher column density gas. 

3.2. CCH column densities 

In the PRISMAS observations we have detected 5 tran- 
sitions of CCH at frequencies around 87.3 GHz and 
524.0 GHz. Wc have first attempted to derive a self- 
consistent LTE model which explains the intensities of 
all the CCH transit io ns using a combinat ion of XCLASS 
( Schilke et all [l999t iComito et all l2005h and MAGIX 
(|Bernst et al.1 l201l[ ). XCLASS generates synthetic spec- 
tra consisting of multiple velocity components as well as 
hyperfinc structure components for all transitions of dif- 
ferent chemical species. XCLASS accesses the molecular 
databases CDMS and JPL and models each molecule with 
the following free parameters: source size, temperature, col- 
umn density, line width and velocity offset relative to the 
systemic velocity of the source, and derives the column den- 
sity corresponding to the different velocity components de- 
tected in absorption and emission in the observed spectra. 



C3H2 30 31±4 0.7±0.1 3.3±0.1 -2.7 ±0.1 

15 78±6 0.9±0.2 1.5±0.2 -4.5 ± 0.2 



The source size refers to the size of the cloud contributing 
to the emission or absorption spectra. The synthetic spec- 
tra generated with XCLASS were fitted to the observed 
spectra using MAGIX, an iterating engine that allows au- 
tomatic minimization to constrain model parameters. We 
have assumed that the hypcrfine lines have the same excita- 
tion temperatures. We note that the fitting procedure can 
involve considerable degeneracy between the source size, 
the excitation temperature and the column density. In the 
absence of more accurate information we have fixed the 
source sizes for MM1 and MM2 to 15" and the size for 
the envelope to 30". Using XCLASS/MAGIX we could not 
obtain an LTE model which consistently explained the in- 
tensities of all the 5 lines of CCH. We found that the 87 GHz 
lines have intensities much higher than the expected values 
based on the 524 GHz intensities under the assumption of 
LTE. We have thus fitted two different LTE models to the 
87 GHz and 524 GHz spectra of CCH. Figure Ega) and (b) 
show the observed CCH spectra along with the LTE models 
fitted using XCLASS and the fit parameters are presented 
in Table [5] The errorbars for the parameters in Table [5] 
represent the range of values over which the best fit model 
does not change significantly. 
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Fig. 4. Spectra of (a) CCH at 87.316 GHz, (b) CCH at 
523.9 GHz, and (c) C 3 H 2 at 85.338 GHz fitted with LTE 
models using XCLASS. (a) also shows the main hyperfine 
component of 87.316 GHz as inset to demonstrate the ve- 
locity structure better. 



From the intensities of the 87.316 GHz CCH lines we 
estimate the LTE temperatures of the —0.5, —2.9 and 
—4.2 kms -1 components contributing to the CCH emis- 
sion to be 110 K, 9 K, and 111 K respectively. The CCH 
transition at 524 GHz shows evidence of only two velocity 
components at —2.8 kms -1 and —4.6 kms -1 , with temper- 
atures of 33 K and 101 K respectively. 

While the temperatures of the envelope derived from the 
LTE models of the 87 GHz and the 524 GHz transitions are 
not too discrepant, the total CCH column densities of the 
envelope (u C cn = —2.7 kms -1 ) estimated from the two tran- 
sitions are 2.3 xlO 15 cm -2 and 2.9 xlO 14 cm -2 respectively, 
hence differing by a factor ~ 8. This discrepancy in column 
densities indicates two possibilities: First, it is likely that 
the excitation of CCH is non-LTE so that, the 87 GHz tran- 
sition with a low critical density is thcrmalized and traces 
the total column density better, but a commensurate num- 
ber of molecules in the higher energy levels probed by the 
524 GHz lines are not present. Secondly, owing to the dif- 
ference in the pointing centers of the beams the 87 GHz 
data could also be tracing the emission from the extended 
ridge in DR21(OH), whereas the higher- J lines at 524 GHz 



primarily trace the DR21(OH) hot core and the envelope. 
Though difficult to quantify, the difference in coupling of 
the emission with the beams at 87 and 524 GHz, which are 
of different sizes and differently centered could contribute 
to part of the discrepancy. For consistent comparison of 
the abundances of CCH, C-C3H2 and C3 abundances in 
the envelope, we use the column densities derived from the 
87 GHz CCH transitions as an upper limit for chemical 
models discussed later. 

3.3. LTE Modelling of c-C 3 H 2 

We have also derived an LTE model for the single 0C3H2 
line using a two-component Gaussian profile correspond- 
ing to the velocities around —2.7 kms -1 and —4.5 kms -1 . 
Fig. 0Jc) shows the observed 0C3H2 spectrum fitted with 
a two-component model using XCLASS. 

3.4. Temperature structure in DR21(OH) 

Based on the results of synthesis imaging of thermal NH3 
and detection of the (7,7) transition of NH3, which is 535 K 
above the ground state, the kinetic temperature of the gas 
associated with MM1 (f lsr = ~4-l kms -1 ) is estimated t o 
be > 80 K (iMangum et"aHll992tlMauersberger et al.lll986l ). 
IWilson fc Mauersbergerl (|1990D derived a kinetic tempera- 
ture of 34 K for this region using NH3 observations at a 
resolution of 40". 

Wc note that the 43" beamsizc for the 524 GHz HIFI 
observations and the 29" beamsize of IRAM 30 m (Fig. 0) 
include MM1, MM2 and the envelope. Hence the observed 
CCH emission is primarily due to the molecular gas sur- 
rounding all these sources which is traced in CO and C 18 0. 
Using CCH and 0C3H2 we derive the temperature of the 
component corresponding to the envelope to be 10-30 K 
and the temperature of the component associated with 
MM1 to be 100-110 K. The data on carbon chain molecules 
are consistent with the presence of high temperatures in 
MM1 and MM2 and agree with the previous determina- 
tions. 

4. Modeling the Observed Chemical Abundances in 
the Envelope 

4.1. The OSU Gas Grain Models 

We model the observed abundances of the three chemical 
species C3, C-C3H2 and CCH, (Table [5]) in the warm enve- 
lope around DR21(OH) using the Ohio State University 
(OSU) gas-grain code with a warm-up (jHasegawa et al.l 
[l99llGarrod fc Herbstll200l . Although slightly shifted in 
velocity from the envelope at —4.2 kms -1 detected in C3, 
for CCH-524 and c-C 3 H2 wc consider the —4.5 kms -1 ve- 
locity component to correspond to the envelope. Since the 
models are primarily for dense star forming regions and the 
low-frequency CCH transition (87 GHz) likely trace an ad- 
ditional extended component, we use the CCH abundance 
derived from the IRAM data as a n upper limit to com pare 
with the outcome of the model. iMotte et al.l (|2007T ) esti- 
mated the average volume density of the entire region to 
be 10 6 cm -3 , while the MM1 and MM2 dense cores reach 
densities larger than 10 7 cm -3 (Table [l}. Thus for the en- 
velope we presently consider models with constant gas den- 
sities n H2 = 5 x 10 5 , 1 x 10 6 , and 5 x 10 6 cm -3 . 
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Table 6. Comparison of observed abundances of the chem- 
ical species in the envelope of DR21(OH) and other sources. 
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Table 7. Parameters for the models explored for C3 chem- 
istry in DR21(OH) 



Model 


T max (K) 


71h 2 (cm" J ) 


1 


30 


5 x 10 s 


2 


30 


1 x 10 6 


3 


30 


5 x 10 6 


4 


50 


5 x 10 5 


5 


50 


1 x 10 6 


6 


50 


5 x 10 6 



Table 8. Parameters assumed to remain constant in 
present models. 



" Estimated u sing N (H 2 ) = 3.5xl0 23 cm" 2 (Table [TJ. 
b ISakai et al.1 [|2009alltft 

c ISakai et al.1 (|2009af i; ICazaux et all (|2003l ); iBisschop et"all 
(120071 s ! 

d Observations not available, Hasscl c t al.1 l|201lT ) model 
(n H2 = 5 x 10 5 cm" 3 ) predicts X(c-C 3 H 2 )« 10" 9 in luke- 
warm corinos (T~ 30 K) and X(c-C 3 H2)~ 10~ 7 in hot cores 
(T~ 100 K). 



Parameter Value Unit 

b Cosmic Ray Ionization Rate £ 1.3 x 10~ 17 s _1 

b Visual Extinction Ay 10 

a,b Grain Radius a d 1 x 10~ 5 cm 

a ' Non-thermal Desorption Yield Orrk 0.01 

a ' b Initial Temperature 71 10 K 

a ' b Warm-up Timescale t^ b 2 x 10 5 yr 

b Warm-up Offset Time t oS 1 x 10 5 yr 



The gas-grain network considers 7166 reactions involv- 
ing a total of 668 gaseous and surface species, where the 
surface species are identified with (s). The physical param- 
eters and i nitial chemical abund ance s adopted are the same 
as us ed by iGarrod et al.l (|2007|) and lHassel et al.l (Tables 1 
& 2; l2008f ). In this approach, a one-point parcel of mate- 
rial undergoes an initially cold period of To = 10 K with 
a duration of t = 10 5 yr, followed by a gradual temper- 
ature increase. A heating timescale of = 0.2 Myr is 
adopted following IGarrod &: Herbstl (|2006t ). where we use 
values for the maximum temperature T max of 30 K and 50 
K, so that the model reaches T max by t sa 0.3 Myr. A to- 
tal of six combinations of T max and hh 2 are considered, as 
outlined in Table [7j In contrast to more complex hydro- 
dynamic approaches (|Aikawa et al.l 120081 ) , here the other 
physical parameters (Table [5]) are assumed to be homoge- 
neous at any given time and from one model to another. 
The advantage of this approach is that it allows a more 
detailed look at the chemical processes and the roles of in- 
dividual reactions. A similar procedure was previously used 
to simulate the formation of hydrocarbon chai ns in the en- 
velope s of low mass protostell ar envelopes by lHassel et al.l 
(|2008h and lHassel et all (|201ll) . where the latter model also 
considered the linear C 3 detected toward W31C and W49N 
(|Mookeriea et al.ll2010D 

The models include photodesorption rates for CO(s), 
N2(s), H2 0(s), and CC>2(s) based on the measurements of 
lOberg et all (120071 . 2009a, b). The adopted processes and 
rates are species-dependent, rather than dependent on a 
single parameter. 

The production of hydrocarbon chains during the warm- 
up has been attributed to a process beginning with 
grain surface chemistry. The elevated temperature allows 
methane to evaporate from the ice mantles surrounding 
dust particles and act as a precursor f or a car bon-c hain rich 
ion-molecule chemistry. lHassel et al.l (Fig. 6; l201lh demon- 
strated that the abundance of C3 is subject to such a pro- 
cess. In particular, the C3 is formed via gas-phase chemistry 



Q IGarrod fc Herbstl (|2006l l: b lHassel et al.l (2008) 



involving atomic C at early times, and depletes sharply at 
t ~ 10 5 yr, after which time most available atomic C is in- 
corporated into CO. In the absence of grain surface chem- 
istry, the C3 will continue to deplete to abundances below 
10~ 12 by 1 Myr. In models including grain surface chem- 
istry, the abundance of C3 is rc-gcncratcd following the re- 
lease of CH4 (s) and subsequent gas phase chemistry during 
the warm-up. The abundances and reactions discussed in 
this section refer to the gas phase. 

4.2. Model results 

The time evolution of fractional abundances is presented 
for C3, CCH and C-C3H2 in Fig. [5l along with the corre- 
sponding observed values. In all the models the warm-up 
occurs between 0.1 and 0.3 Myr starting from an initial tem- 
perature o f 10 K, simil ar to the te mper ature profiles that 
appear m lHassel et all (Figs.2&3; I2008D . We observe the 
following general trends in the abundance profiles shown in 
Fig. [5] For the lower temperature models (1-3), the molec- 
ular abundances attain values comparable to the observed 
abundances for a much longer period of time as compared 
to the higher temperature models (4-6). Further, for both 
C3 and C-C3H2 in Models 4-6, the maximum abundance 
decreases with increasing values of ?ih 2 • 

Fig [5] shows that all six models produce CCH to a peak 
abundance of A = 0.2 — lxl0~ 8 during the warm-up. The 
abundance depletes more rapidly for Models 4-6, while the 
abundance remains elevated for a longer period for Models 
1-3. Moreover higher densities produce a slightly smaller 
CCH abundance following the warm-up. The C3 abundance 
depletes more rapidly for Models 4-6, and smaller abun- 
dances are produced for larger density during this period. 
All the models produce a similar abundance of 0C3H2 in 
the range of X = 0.5 — 1 x 10~ 8 , following the warm-up. 
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Fig. 5. Comparison of observed abundances in the envelope (shown by dashed lines) with the results of chemical models 
involving warming up of grains for different combinations of hydrogen density (?t-h 2 ) an d final warm-up temperatures. 
The model calculates abundances relative to hydrogen, and the observed column densities of the different species have 
been converted to fractional abundances assuming N(ti2) = 3.5x 10 23 cm~ 2 . In all the models the warm-up occurs 
between 0.1 and 0.3 Myr starti ng from an initial temperature of 10 K, similar to the temperature profiles that appear in 
lHassel et al.1 (Figs.2 & 3: l2008l) . 



The temperature and density effects seen in the previous 
two species also appear for C-C3H2 but to a lesser extent. 

The primary reason for the difference in abundance of 
the chemical species considered here, for T max = 30 K and 
50 K models is that the higher temperature will release ad- 
ditional species from the icy grain mantles to the gas. This 
leads to additional, competitive reaction pathways with the 
overall effect of a comparative depletion of abundance of 
these particular species. A similar effect is visible at the 
onset of the warm-up, when the abundances of the three 
species deplete sharply just prior to the steep increase at- 
tributed to the warm-up. 

4.3. Major Chemical Reactions 

Table summarizes the most important chemical reactions 
occurring in the gas-phase during warm-up as determine d 
by the gas-grain chemistry model (jGarrod fc Herbstll2006l ) . 
We emphasize that the species produced here also result 
from other ion-molecule reactions, and that Tabic |H] lists 
only the most direct and salient mechanisms. The liber- 
ated CH4 reacts with C + to form the intermediate species 
C2HJ an d C2H3" via reactions (1) and (3). Subsequently 
reaction of H2 with C2HJ in (2) produces C2H4 ■ The dis- 
sociative recombination reaction of an electron with these 
intermediate ions then produces either C2H2 by reactions 
(2) or (5). Further reaction of CH4 with the intermediate 
ions can produce CsH^ via reactions (6) and (7), while the 
dissociative recombination of CsH^" produces C3H3 via re- 
action (8). Reactions of C and C + with C2H2 (9,11) produce 
C3H and C3H + . Further reaction of these products (12,13) 
produce C3H2 , and the reaction of CH4 with C3H2 (14) 
produces C3HJ. The C + and C result from partial des- 
orption of CO(s) and subsequent reaction with Hc + during 
the warm-up. These represent the major precursors for the 
observed species, and provide the basis to examine each 
species individually. 

For the chemical species considered here we identify the 
primary formation reactions as follows: CCH is primarily 
formed by the dissociative recombination of C2H^" in re- 
action (17). The evolution of C3 proceeds in a somewhat 
different manner. Ion-molecule chemistry forms a relatively 



large abundance of this species at early times, followed by 
a depletion to grain surfaces. This depletion happens at 
slightly earlier times for models with larger values of den- 
sity. Although C3 is frozen on the dust grains, during the 
warm-up, C3 is not liberated directly from the grains, but 
rather is produced in a reaction sequence starting with CH4, 
whilst the Cs(s) is more likely to participate in surface re- 
actions than to desorb. The reaction of C + C2H2 pro- 
duces C3 through reaction (19). This mechanism and the 
secondary pathways result in the formation of a peak in 
the abundance during and after the warm-up period. The 
formation of C-C3H2 can proceed via three pathways. The 
first mechanism is reaction (21), in which C3H3 reacts with 
H. The dissociative recombination of C3HJ (22) is a second 
mechanism. Alternatively, reaction (23) of CH with C2H2 
is another important formation pathway. 

4.4. Best fit chemical model 

In order to compare model results with observations on a 
more quantitative basis, w e compute the mean confidence 
level (jGarrod et all l2007t) during the warm- up period. In 
this q uasi-statistical analysis introduced by iGarrod et al.l 
(|2007t ). the confidence level, Ki, for agreement between com- 
puted and observed abundances is defined as 



erfc 



|log(X i )-log(X o6a , i )| 



(1) 



The index i refers to individual chemical species un- 
der consideration. The mean of the individual values of «j 
thus computed is denoted by n avg . We consider the com- 
puted abundance to "fit" the observation if the abundance 
agrees within an order of magnitude or better, and thus 
assume a = 1 so that one standard deviation corresponds 
to agreement of one order of magnitude. The related confi- 
dence level for such a factor of 10 agreement is Ki — 0.317, 
and the value increases up to Ki = 1 for closer agreement. 
The accuracy of this analysis and the robustness of K avg 
against individual values of Ki improves with increasing size 
of the considered dataset. Since we have considered only 
three chemical species, we use the ft-based analysis primar- 
ily as a guideline to identify models which best reproduce 
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Table 9. Dominant reactions of observed and intermediate 
species following CH 4 sublimation. The observed species 
arc identified in bold type. 
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Fig. 6. Mean confidence level K avg cal culated for the chem - 
ical models as explained in the text (jGarrod et al.l l2007h . 



the data. The best fit corresponds to the maximum average 
confidence level. 

Fig. [5] shows the K avg values for each model that we 
have considered here, along with horizontal lines indicating 
agreement within a factor (/) of 2, 3, 4, and 5. The primary 
interpretation of this diagram is that all six models reason- 
ably fit the observed abundances during the period follow- 
ing the warm-up. Models 1-3 fit the observation set within 
a factor of 2-4 or better for the time period of t ~ 0.7 — 3 
Myr for all three density values. Models 4-6 fit the obser- 



vation set within a factor of 3-5 or better for the period of 
t ~ 0.6 — 0.8 Myr for Models 4 and 5 and for the earlier 
period t ~ 0.3 - 0.7 Myr for Model 6. Model 3 gives the 
best agreement with observations in terms of the value and 
duration of the K avg function, while Model 2 is the second- 
best at meeting these criteria. This analysis quantitatively 
asserts the level of model agreement described in Fig. [5l 
and broadly shows that the best sustained agreement of 
the models occurs in the time following the warm-up rather 
than earlier or later. 



5. Discussion 

5.1. Lower gas phase abundance of C3 in MM1 

The angular resolution of the spectroscopic observations 
used in this paper lies between 12 and 40", so that MM1 
and MM2 are not resolved spatially. However, due to their 
markedly different velocities, MM1, MM2 and the envelope 
are easily discernible in the spectra. Interestingly, all C3 
spectra show primarily two velocity components associated 
with MM2 and the envelope, with absorption due to MM1 
being detected only in Q(4). In contrast, the chemically 
related species CCH and C-C3H2 that we have considered 
here show components corresponding to only MM1 and the 
envelope. Thus there is an intrinsic difference in the dis- 
tribution of C3 compared to the distribution of the other 
species. The abundance of C 3 in MM2 (1.2xl0~ 9 ) is only 
twice that in the envelope (6.3xl0 -10 ). Based on Tabled] 
the total column density of MM1 is ~ 60% of the column 
density of MM2. This suggests that if the relative abun- 
dances of C3 in MM1 and MM2 were similar, the total col- 
umn density of C3 in MM1 would have been 2.5 times the 
column density in the envelope so that the absorption dips 
due to MM1 would be similar to or stronger than the dip 
due to the envelope. It may be argued that the continuum 
backgrounds for the two components MM1 and the enve- 
lope may be different. However the continuum from MM1 
is expected to be stronger than the continuum from the 
envelope. All of the above arguments indicate that C3 is 
preferentially depleted or destroyed in MM1. 

One major difference between MM1 and MM2 lies in 
the fact that MM1 has already shown evidence for at least 
one high temperature hot core, whereas MM2, though more 
massive shows no such components. The chemical network 
i ncludes among others the rea ction C3 + H 2 — > C3H + H 
(|Pineau des Forets et al.|[l988|) . which becomes a primary 
destruction pathway for C3 when the temperature exceeds 
« 80 K. As a result, models predict that C3 will be de- 
pleted in favor of 0C3H2 in hot core conditions, as C3H is 
also destroyed by reaction with H2 . Additionally, the chem- 
ical models presented here also show that the C3 abundance 
is reduced at larger densities. MM1 is definitely at a higher 
density than the envelope (Table [I} , so that the reduced 
abundance of C3 in MM1 relative to MM2 could also indi- 
cate that MM2 has a density lower than MM1. We note that 
this contradicts the densities derived from dust column den- 
sities (Table [T|) , which is not completely unexpected since 
the densities derived from the column densities are beam 
averaged and need not always reflect the local densities in 
a region. 

However, an alternate possibility arises from the fact 
that the source intrinsic continuum opacity can re-fill 
the absorption to a significant extent, so that the col- 
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umn densities derived here are all primarily lower limits 
(jMookeriea et al.ll2010t) . It can thus be envisaged that in 
the case of MM1, which has a stronger continuum emis- 
sion, it is a combination of the geometry of the source and 
the re-filling of the absorption dip which reduces the "visi- 
bility" of C3. A proper evaluation of the effect of the con- 
tinuum opacity on the absorption depth is possible only by 
constructing a complete radiative transfer model of the en- 
tire region with accurate temperature and density profiles, 
which are not yet available. 

Finally, since MM2 is estimated to be younger than 
MM1 and c-C 3 H 2 , and CCH emission from MM2 are not 
detected, the possibility that the observed C3 abundance in 
MM2 is due to the pre-warm-up gas-phase chemistry can 
also not be ruled out. 

5.2. Comparison with C3 in diffuse clouds 

C3 has been found to have an al most one-to-one cor - 
relation with C2 in diffuse clouds (jRoueff et al.l l2002t) . 
Further, the abund ance of C3 in diffuse clouds, (3-6) xlO~ 9 
(|Roueff et al.ll200"2h , is ten times larger than the abundance 
of C3 deriv ed in the d ense star forming environment of 
DR21(OH). lOka et al.l (|2003h had explained the observed 
correlation between C2 and C 3 in diffuse clouds, in terms of 
a direct pathway of formation of C3 from C2 (see their Fig. 
4). In the previous section, we examined the formation of 
C3 during a warm-up period, and noted a relatively larger 
abundance of this species prior to t = 10 5 yr, the so-called 
"early time" . We now explore the chemistry of C3 in dif- 
fuse regions and during the early time period in dense gas 
regions, to examine the differences (if any). 

The first step of the transformation of C2 to C3 in diffus e 
clouds is the photoionization of C2 to (|Oka et al.ll2003t) . 
In contrast, the models for the dense star forming regions 
assume darker (Ay = 10) conditions, so photoionization is 
not an effective process. 

In dense models, the process begins with the radiative 
association of C + + H2 — > CHj . Next, CH^ reacts with 
H2 to form CHj. This is followed by the dissociative re- 
combination of CH3" and CH^ as competitive reactions to 
form CH, and then an ion- n eutral reaction w ith C + to form 
C2 • As shown in Fig. 4 of (|Oka et al.ll2003l ). the pathway 
then follows along the upper (C2H) branch to form C3. 
The C3 is destroyed by processes that form C4, C4 , and 
C5, which react further to primarily form C3 in a cyclic 
process. The abundance of C3 grows until C + is depleted 
at t ~ 10 TOO yr, depending on the density, after which 
it reaches a plateau until the C3 is accreted onto grain 
surfaces. Dissociative recombination reactions with C2H" 1 ", 
C2H2 ' an d along the pathway form C2. The C2 then re- 
acts with O to form CO, and as a result, the abundance falls 
to X(C 2 ) ~ IP" 9 , far below the X(C 2 ) : X(C 3 ) = 40 : 1 
measured bv lOka et al.l (|2003h for diffuse clouds. 

Thus, C2 is not the starting point for C3 formation in 
dense cloud models at early times, and also does not un- 
dergo a cycling process like that of C3. As a result, the cor- 
relation between C2 and C3 abundance is not predicted by 
the dark cloud models. Rather the C2 abundance evolution 
more closely follows the evolution of CCH from early times 
until the warm-up begins than that of C3. This comparison 
indicates that the chemistry of C3 formation in the enve- 
lope of DR2I(OH) is definitely not the same as in diffuse 



clouds. Further, the early time abundance of C3 for dense 
regions exceeds the observed abundance by three orders 
of magnitude, effectively excluding this formation route as 
well. 



5.3. Warm Carbon Chain Chemistry in the envelope 

Comparison of the present model results and observations 
indicates that the observed composition can be simulated 
with ion-molecule chemistry following a moderate warm- 
up to 30-50 K. This also means that C3 formation in star 
forming regions can be explained in terms of the moderate 
temperature (30 K) gas-phase chemical reactions starting 
from CH4 evaporated from the grain mantle only, and does 
not require photochemistry of PAH molecules. 

The observed abundances of C3 are consistent with for- 
mation during a warm -up model, or "W arm Carbon Chain 
Chemistry" (WCCC) (ISakai et al.ll2f)08h . and related to the 
abundances of CCH, C-C3H2. This suggests that the enve- 
lope conditions are more like those of "lukewarm corinos" 
surrounding low mass protostars than those of hot cores. 
In other words the C3 is not residing in a hot core. The 
abundances of the additional species chemically related to 
C3, which we consider for the case of DR2I(OH) provided 
further insight regarding the envelope conditions. 

Table [6] summarizes the observed abundances of the dif- 
ferent species in DR2f(OH) as well as those available in 
the literature for two types of sources: the lukewarm cori- 
nos L1527 and B228 and the hot corino IRAS 16293-2422. 
We find that while C3 has not been observed in any of the 
other sources, the abundances of the other three species in 
DR21(OH) are closer to the abundances observed in luke- 
warm corinos. Further, the chemical model that best rep- 
resents the observed abundances corresponds to a T max of 
30-50 K, a temperature that is much lower than the values 
expected in hot cores/corinos. This suggests that the enve- 
lope of DR21(OH) is chemically closer to a WCCC region 
such as a lukewarm corino envelope than to a hot core. 
Further, the pre-depletion C3 abundance (at t < 10 4 yr) 
exceeds the observed values by a large factor. This implies 
that the C3 abundance is primarily maintained via ion- 
molecule chemistry in the gas phase after CH4 is desorbed 
from the surface of the dust grains, and does not have a 
one-to-one correlation with C2 as found in diffuse clouds. 

5.4. Chemical & Dynamical Ages of the Region 

The chemistry of the envelope surrounding the embedded 
sources MM1 and MM2 in DR21(OH) appears to have oc- 
curred over a period between (0.7-3) Myr including a cold 
period, warm-up, and extended time at T = 30 or 50 K. 
This time period appears to be on the longer side of the 
evolutionar y timescales of massiv e star forming cores. As 
detailed by ICsengeri et al.l (|201lf ) the life time of the gas 
in massive dense cores is determined by the crossing times 
or the local free-fall times for the molecular dense cores 
which is (5-7) x 10 yr. However the large scale flows asso- 
ciated with Cygnus X (and DR21 in particular) are massive 
enough to continuously replenish the mass of these cores al- 
lowing them to remain "active" for a much longer period 
of time and the dynamical ti me of the most massiv e sub- 
filament in DR21 is ~ 2 Myr (jSchneider et al.ll2010D . Thus 
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the chemical age of ~ 1 Myr is consistent with the dynam- 
ical age of the region. 

6. Summary 

In the DR21(OH) region, using the high velocity resolution 
of HIFI, we detect absorption due to C3 from the envelope 
and the less evolved core (MM2), and a rather weak signal 
associated to the most massive and hottest region (MM1). 
The abundance of C3 in the envelope of the hot core associ- 
ated with DR21(OH) has been consistently explained along 
with those of other species formed in the same chemical 
network using a chemical model involving the warm-up of 
grains. The formation mechanism of C3 required to explain 
the observed abundances in dense star forming regions is 
the gas-phase reaction of CH4 desorbed from the grain sur- 
faces that have been warmed up. It is not formed in a cycle 
with C2 as the starting point nor is it produced by the pho- 
todestruction cascades of PAH molecules. For the envelope 
of DR21(OH) the chemical models which best explain the 
abundances of C3 and other chemically related species cor- 
respond to Models 2 and 3, nn 2 =l-5xl0 6 cm -3 for T max 
= 30 K , for a time period between 0.7-3 Myr and Model 6, 
iih 2 =5 x 10 6 cm~ 3 for T max = 50 K, for a time period be- 
tween 0.6-0.8 Myr. The upper limit of the timescale, though 
much larger than the typical lifetimes of hot (warm) cores, 
is not unrealistic in the case of DR21(OH) where the dy- 
namic age of the most massive filament has been estimated 
to be 2 Myr. 
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